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CN ; Abstract 

This paper is concerned with the 1-D compressible Euler-Poisson equa- 
tions with moving physical vacuum boundary condition. It is usually used 
to describe the motion of a self-gravitating inviscid gaseous star. The lo- 
cal well-posedness of classical solutions is established in the case of the 
, adiabatic index 1 < 7 < 3. 



1 Introduction 



< 

f~i ' The motion of self-gravitating inviscid gaseous stars in the universe can be de- 

■ scribed by following free boundary problem for the compressible Euler equations 
coupled vifith Poisson equation: 

+ V • (pu) =0 in n{t), (1.1) 

p[ut + u-Vu]+VP ^ pV(f> in fl{t), (1.2) 

: -A0 = 4^P5 on n{t), (1.3) 

c^: Hm) = u, (1.4) 

Oi ip,u) = {pa,uo) on n{0). (1.5) 

■ The open, bounded domain Q{t) C denotes the changing domain occupied by 
^—^ ' the gas. r(i) :— dQ{t) denotes the moving vacuum boundary, i'{T{t)) denotes 

the velocity of T{t). The density p > in n{t) and p = in M"^ \ il{t). u denotes 
the Eulerian velocity field, p denotes the pressure function, and p denotes the 
density of the gas. cf) is the potential function of the self- gravitational force, and 
I g is the gravitational constant. The equation of state for a polytropic gas is 

^ ■ given by: 

P = C^pT for 7 > 1, (1.6) 

where C-y is the adiabatic constant which we set to be one. For more details of 
the related background, see for instance, [1]. 

With the sound speed given by := y/dP/dp, and with cq = c(-,0), the 
condition 

-00 < ^ < on r (1.7) 
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defines a "physical vacuum" boundary, wliere N denoting the outward unit 
normal to the initial boundary F := 917(0). This definition of physical vacuum 
was motivated by the case of Euler equations with damping studied in [5] [5] and 
the physical vacuum behavior can be realized by some self-similar solutions and 
stationary solutions for different physical systems such as Euler equations with 
damping. For more details and the physical background of this concept, please 
see [S][H][Tn]- 

The local existence theory of classical solutions featuring the physical vac- 
uum boundary even for one-dimensional compressible Euler equations was only 
established recently. This is because if the physical vacuum boundary condition 
is assumed, Euler equations becomes a degenerate and characteristic hyper- 
bolic system and the classical theory of hyperbolic systems can not be directly 
applied. In [T^], Jang and Masmoudi consider the one-dimensional Euler equa- 
tions in mass Lagrangian coordinates. Local existence was proved using a new 
structure lying upon the physical vacuum in the framework of free boundary 
problems. Independently of this work, in [2], Coutand and ShkoUer constructed 
-ff^-type solutions with moving boundary in Lagrangian coordinates based on 
Hardy inequalities and degenerate parabolic regularization. 

In this paper, we will focus on the 1-D case for the system (|1.1[) - (|1.5|) with 
the physical vacuum condition: 



Pt + {Pu)r, 


= 


in 


m. 


(1.8) 


p[ut + UUjj] + {p''),, 


= P(l>v 


in 


m, 


(1.9) 


't'riri 


^Cp 


on 


m. 


(1.10) 


{P,u) 


= (po,wo) 


on 


m, 


(1.11) 




= u, 






(1.12) 


0< f 
or] 


< -l-oo 


on 


r, 


(1.13) 



where /(O) — I = {Q < r] <1} and F := dl and prove the local existence result 
for it (The 3D case will appear soon). 

Our main result is the following theorem: 

Theorem 1.1. (Local wellposedness) For 1 < 7 < 3, assume that initial 
data po > in I, Mq < 00 (defined in i2.26]) ). and the physical vacuum condi- 
tion il.lc)]) holds. Then there exists a unique solution to i2.3\) - ![^7E\) (and hence 
il.8\) ~ n.l2\) ) on [0,r] for some sufficiently small T > such that 

sup E{t) < 2Mo. (1.14) 

te[o,T] 

The local wellposedness result for the corresponding Euler equation was 
proved by Coutand and ShkoUer in [2]. Motivated by [2], we also use the La- 
grangian coordinates to reduce the original system to that in a fixed domain. 
In our problem, we have the extra potential force term 0,, in (|1.9p . To handle 
this term, we will give an explicit formula for it and show that it is a function 
of pq. Then we construct the approximate solution in two steps. Firstly, using 
Galerkin scheme to find the solution to linearized problem of the degenerate 
parabolic regularization. In this process, we would make fundamental use of 
the higher-order Hardy-type inequality introduced in [5] (we would give a new 
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proof). But we will define an intermediate variable which is different to the 
one used in [2]. By using our intermediate variable, the improvement of the 
space regularity for the solution of linear problem will be easy and clear with 
less computation. Secondly, using fixed-point scheme to get the solution to the 
degenerate parabolic regularization. Last we would derive a priori estimates for 
the approximate solution. This part is more or less similar to that in P]. For a 
self-contained presentation, we will still carry out the proof in Section [51 Below 
we will mainly focus on the case of 7 = 2. The general case for 1 < 7 < 3 is 
treated in Section [51 

Now we briefly review some related theories and results from various aspects. 
For Euler-Poisson equations, the existence theory for the stationary solutions 
has been proved by Deng, Liu, Yang, and Yao in [TT]. For Navier-Stokes- 
Poisson equations, Li, Matsumura and Zhang studied optimal decay rate for 
the system and Zhang, Fang studied global behavior for spherically symmetric 
case with degenerate viscosity coefficients in ^B] . 

For compressible fluids, Makino proved the local-in-time existence of solution 
in [18] with boundary condition p = for some non-physical restrictions on 
the initial data. And Lindblad proved the local-in-time existence with vacuum 
boundary condition P = for general case of initial data with the main tool 
which is the passage to the Lagrangian coordinates for reducing the original 
problem to that in a fixed domain in 17 . And H.L. Li, J. Li, Xin [23], Luo, 
Xin, Yang [53], Xin[5S] also did many works on compressible Navier-Stokes 
equation with vacuum. 

For incompressible flows, Wu solved local well-posedness for the irrotational 
problem, with no surface tension in all dimensions in [14) and |15) . Lindblad 
proved local existence of solutions general problem with no surface tension, 
assuming the Rayleigh- Taylor sign condition for rotational flows in |16| . For 
the problem with surface tension, B.Schweizer proved existence for the general 
three-dimensional irrotational problem in [19) . And we also mention the works 
by Ambrose and Masmoudi [20], Coutand and Shkoller 4 , and P. Zhang and 
Z.Zhang [H]- 

This paper is organized as follows: In Section [2] we formulate the problem 
in Lagrangian coordinates. In Section [3] we present some lemmas will be used. 
In Section [3]-Section[Sl we introduce a degenerate parabolic approximation and 
solve it by a fixed-point method. In Section[n]-[71 we derive the a priori estimates 
and prove the local well-posedness for 7 = 2. In Section[Sl we discuss the general 
case for 1 < 7 < 3. 

2 Lagrangian Formulation 

Here, we denote rj as Eulerian coordinates and denote x as Lagrange coordinates, 
which means r](x, t) denotes the "position" of the gas particle x at time t: 



where o denotes the composition [uo r]]{x,t) — u{ri{x,t),t). We also have: 



dti] = u o rj for i > and 77(2;, 0) = x, 



(2.1) 



V = uo rj (Lagrangian velocity) , 

f — p° rj (Lagrangian density), 

<^ = (j)orj (Lagrangian potential field) . 



(2.2) 
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2.1 Fixing the domain and the Lagrangian version of the 
system 

Noticing (|2.ip and (I2.2p . the Lagrangian version of system (I1.8p ~ (ll.l2[) can be 
written on the fixed reference domain / as 

ft + f^v = in /X (0,r], (2.3) 

fvt + -^if) = fAr^ in / X (0,r], (2.4) 

-{-^)''i' = Cf in /X (0,r], (2.5) 

{f,v,v) ^ {Po,uo,e) in/x{i = 0}, (2.6) 

where e{x) = x denotes the identity map on /. 
By conservation law of mass, we have 

f = poT] = po/r]^. (2.7) 

Hence, the initial density function po can be viewed as a parameter in the Euler 
equations. 

Since po > in /, (jl.l3p implies that for some positive constant C and x G / 
near the vacuum boundary F, 

po > Cdist(x,r). (2.8) 
Hence, for every x €i I, we have: 
dpo 



> C when d{x,dl) < a, (2.9) 
Po>Ca>0 when d{x, dl) > a (2.10) 



for some a > and a constant Ca depending on a. 

In summary, we write the compressible Euler-Poisson System as 

PoVt + {pl/ril)x^ Po^x/Vx in/x(0,T], (2.11) 

-i<^x/Vx)x=Cpo in/x(0,T], (2.12) 

{v,r])^{uo,e) in / x {t = 0}, (2.13) 

po=0 on r, (2.14) 

with pq > Cdist(a;, F) for x (z I near F. 

2.2 The formula for potential force 

Now wc try to give an explicit formula for the potential force (f>r] in (II. 8p and 
corresponding term ^^/Vx in (|2.11l) . Set I{t) = {a{t),b{t)) = (77(0, i), 77(1, i)). 
First, for every i, we can directly get 

Mv,t) = - f\{y,t)dy + Mit). (2.15) 

Ja{t} 
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It is reasonable to assume that oo)| ~ |(/),,(+oo)|. Since the gas only 

occupied bounded interval, the force of gas produced in — oo and -\-oo can be 
regarded as the same large (see [S])- Noticing — —Cp < 0, so we get: 



+ C30 

.^(+(x)) = -^/ p{y,t)dy, (2.16) 



1 
2 

1 
2 



-\-oo 

-ni-^) = ^l p{y,t)dy, (2.17) 



+00 ra{t) \ 

p{y,t)dy- / p{y.t)dy]. (2.18) 
Since pir], t) > in I{t) and p(?7, t) = when 77 < a(t) or 77 > fe(i), we have 

Ja(t) ^ Ja(t) 

Then we transform the formula (|2.19p to Lagrange variables: 

p{y,t)dy + - / piy,t)dy 

a{t) ^ Ja(t) 

ri{x,t) /•»7(l,t) 

p[ri[z, t), t)dri{z, t) + - I p{ri{z, t), t)dr]{z, t) 
))(o,t) ^ Jn(o,t) 

1 

f{z,t)r]^dz + - f{z,t)r]^dz. 

(2.20) 

With (1^ . we can finally get 

F:^^ = - r po{y)dy + \ ( po{y)dy. (2.21) 
Vx Jo ^ Jo 

Remark 2.1. If po E C", t/ien F e C^+", we wiZ^ see that this regularity 
is important for the case 7 7^ 2 m Section 

With formula (|2.2ip . we can write the whole system as 



Povt + {plhDx = PqF 
Po = 



in I X (0,T], (2.22) 
in I X {t = 0}, (2.23) 
on r, (2.24) 



with Po ^ Cdist(a;,r) for x E I near F. 
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2.3 The higher-order energy function 

The higher-order energy function is defined as follows: 

4 2 

+ w^dtdlvit, + Wpldtdlvit, 

+ \\Vp^dfdMt,.)\\lm) + \\pld!dlv{t,.)\\l.^jy 



(2.25) 



Let P denotes a generic polynomial function of its arguments whose meaning 
may change from line to line. Let 

Ma = P{E{0;v)). (2.26) 

3 Weighted Spaces and A higher-order Hardy- 
type inequaUty 

3.1 Embedding of a Weighted Sobolev Space 

Using d to denote the distance function to the boundary F, and letting p = 1 or 
2, the weighted Sobolev space iJ^p(/), with norm given by ( Jj d{x)P{\R{x)\'^ + 

\dxR{x)\'^)dx) ^ for any R g H^p{I), satisfies the following embedding: 

HUI)^H'~Hi). (3.1) 
So that there is a constant C > depending only on / and p such that 

ll^lli-p/2 < C J^d{xn\Rix)\' + \d,Rix)\')dx. (3.2) 

See, for example, Section 8.8 in [7]. 

3.2 A higher-order Hardy-type inequahty 

The following two lemmas can be found in [2^ . We will use Lemma 13.11 to 
construct the approximate solution in Section [5] and use Lemma 13.21 to obtain 
estimates independent of the regularization parameter defined in Section 21 

Lemma 3.1. Let s > 1 be a given integer, and suppose that 

ue H'{I)nH^{I), (3.3) 
and d is the the distance function to dl, we have that ^ G H''^^{I) with 

^ <C\\u\\h^. (3.4) 
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Proof. Let u e H'{I) n H^{I), then for < m < s: 
where 

/ = ^ C;^a™-'=u(-l)'=A;!d'"-'=. (3.6) 

k=0 

With simple calculation, we can get 

d^f = d^+^u{-l)"'m\d'". (3.7) 
Now using the fundamental calculus theorem, when < x < 5, we have: 

fix) = m+ r dJiy)dy (3.8) 
Jo 

= x I d^f{9x)de (3.9) 
Jo 

= (-l)'"m!a;'"+^ / d^+^u{ex)e'^de. (3.10) 
Jo 

Similarly, when ^ < a; < 1, we have: 

f{x) = (-l)'"+im!(l - x)"'+^ [ d^+\{l - 6»(1 - a;))6»"rf6». (3.11) 

Jo 



Then: 



ii5r(^)iii^ = 



- 1 
Jo Jo 



(3.12) 



1 rl 



+ I [(-1)"+^to! / d^+^u{l-0{l-x))d'"de]^dx 
2 

1,2(7) • 







In this way, we finally get: 

u 

W^Wh^ <C\\u\\hs+i. (3.13) 

□ 

Lemma 3.2. Let k > and g e L°°{0,T; H%I)) be given, and let f € 
H'^{0,T;H'{I)) be such that 

f + Kft=g in (0,T)x7. (3.14) 

Then, 

II/I|l~(0,T;//»(/)) < Cmax{ 1 1/(0) 11^.(7), ||5||L=o(0,T;i/=(/)}, (3.15) 

where C is independent of k. 



4 The degenerate parabolic approximation of the 
System 

4.1 Smoothing the initial data 

For the purpose of constructing solutions, we will smooth the initial velocity 
field uq and density field po while preserving the conditions po > in / and 
(EH) firstly. 

For K > 0, let < ai^{x) E C^(M) denote the standard family of moUifiers 
with spt{aii) — {a;||x| < k} and let Ej denote a Sobolev extension operator 
mapping to iJ'*(]R) for s > 0. 

Now we set smoothed initial velocity filed Uq as: 

^^0 = "l/|lnK| *^^/(wo), (4.1) 

and smoothed initial density function Pq is defined as the solution of the elliptic 
equation: 

d^,p^^d^,[a,/\i^^l*Ei{po)] in /, (4.2) 

= on r. (4.3) 

So for sufficiently small k > 0, Uq,Pq G C°°{I), pg > in /, and vacuum 
condition (|2.8I) is preserved. Details can be found in [2]. 

From now on, we will denote ujj by uq and Pq by po for convenience and it is 
easy to show that Theorem 1 1 . II holds with the optimal regularity by a standard 
density argument. 

4.2 Degenerate paraboHc approximation 

For notational convenience, we will write 

and similarly for other functions. Now for k > 0, we consider the following 
nonlinear degenerate parabolic approximation of the compressible Euler-Poisson 
System ^^TI^ - ^^:^ : 

2 

Povt + = PoF + Kiplv'Y in/x[0,r], (4.5) 

(w,77) = (uo,e) in / x {t = 0}, (4.6) 

po = on r (4.7) 

with po{x) > Cdist(j;,r) for x € I near F. We emphasis that the data (poi^o) 
has been smoothed as in Section (|4.ip . 

We will first obtain the existence of the solution to (|4.5p - (|4.7p on a short 
time interval [0, T^] (with possibly depending on k). Then we will show that 
the time of existence does not depend on k via a priori estimates in Section [5] for 
this sequence of solutions independent of n. Then the existence of a solution to 
the compressible Euler-Poisson system is obtained as the weak limit as k — >■ 
of the sequence of solutions to (I4.5p - (I4.7I) . 
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5 Solving the parabolic k, - problem by a fixed- 
point method 

5.1 Assumption on initial data 

Using the fact that ri{x,0) = x and F — — po{y)dy + Jq Po{y)dy, the 
quantity Vt\t=o for the degenerate parabohc K-problem can be computed using 



Ul 



Vt 



t=0 



Poiy)dy + ^ [ Po{y)dy + —[pIv']' - —{^)'\ 
2 J-oo Po Po V J 



(5.1) 



1 

Po{y)dy+- 



Inductively, for aU fc > 2, fee N: 



Po{y)dy H [Po"o]' - '^P'o 

Po 



k-l 



Po 



■[pi 



Po\"' 



(5.2) 



These formulae make it clear that each d^v\t=() is a function of space-derivates 
of uq and pa. 



5.2 Functional framework for the fixed-point scheme 

For T > 0, we shall denote by Xt the following Hilbert space: 

Xt = {v\v e W^^^{0, T; H\I)) n W^''^{0, T; H^{I)); 
pov e W^'^iO, T; H^{I)) n W^^^{0, T; H^{I))}. 
which is endowed with its natural Hilbert norm: 

Iklli-T =lbllw'5.2(0,T;_ffl(-f)) \\^\\w'^.^(0,T;H^{I)) 

+ IIPow|Ih'5.2(o,T:H2(/)) + ||Pow||w'4,2(o,T;_ff3(7))- 



(5.3) 



(5.4) 



For M > given sufficiently large, we can define the following closed, 
bounded, convex subset of Xt'. 



Ct{M) ^{vCXt: d^v[. 



0,l,2,3,4,5,6,||t;||^^<M}. (5.5) 



which is indeed non-empty if M is large enough which would be determined by 
initial data. Henceforth, we assume that T > is given independently of the 
choice oiv € Ct{M), such that 



ri{x,t) 



v{x, s)ds 



(5.6) 



is injective for t e [0, T], and that ^ < r]'{x, t) < | for t e [0, T] and x el. This 
can be achieved by taking T > sufficiently small: with e{x) — x, notice that 



1^1 = II f v'(x,s)ds\\H^ < VtM. 
Jo 



(5.7) 
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We will apply the fixed-point methodology in AV to prove the existence of a 
solution to the k- regularized parabolic problem (I4.7P . 

Finally, we define a polynomial function Ao of norms of the non-smoothed 
initial data uq and po as follows: 

^fo^P.{\\po\\L-,\\uo\\L-), (5.8) 

where is a generic polynomial with coefficients dependent on powers of | In k|. 

Using the properties of the convolution (|4.ip and (|4.2p . Vs > 1, Vfc G 
1,2,3,4,5,6, the quantities defined in (15. 2p (using the smoothed initial data 
Uq and Pq) satisfies: 

\\uk\\H^ < P{M\\hs+.A\K\\h^^-^) < C.P^QIpoWl^-, ||"o||lO < A/'o- (5.9) 

5.3 Linearizing the degenerate parabolic ^-problem 

For every v € Ct{M), we define tj = x + jQv{x,T)dT and consider the linear 
equation for v: 

povt - k[pIv']' ^ -[||]' + PoF, 

V{-,Q) = Uq, 

where F is defined in (|2.2ip . 

In order to use high-order Hardy-type inequality, it will be convenient to 
introduce the new variable X — pqv, which belongs to Hq{I) (can be seen 
below). Here we choose a different variable X from that used by Coutand 
and Shkoller in [2], which would simplify the process of improving the space 
regularity for solution of (I5.10p . 

By a simple computation, we can see that (jS.lOp is equivalent to 

v,^n^[plv'r = -^[^]' + F, (5.11) 



Po Po rj ■ 



and hence that 



X 



— -kX" + k^X^G in/x[0,T], (5.12) 
Po Po 

X = on r X [0,r], (5.13) 

X\t=o = pqUo in /, (5.14) 



where 



G = F+ ^{^y = - r po{y)dy +]- j p,{y)dy + |,(^)'. (5.15) 



n n Jo 2 7o ?7' 77 

We shall therefore solve the degenerate linear parabolic problem (|5.12p - 
(|5.14p with Dirichlet boundary conditions, which (as we will prove) will surpris- 
ingly admit a solution with arbitrarily high space regularity (depending on the 
regularity of G on the right-hand side of (|5.12p and the initial data of course), 
and not just an _ffQ(T)-type weak solution. After we obtain the solution X, we 
will then easily find our solution v to (|5.10p . 
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In order to apply the fixed-point theorem, we shall obtain estimates for v with 
certain high space-time regularity. Here, we study the sixth time-differentiated 
problem and define the new variable 

Y = dfx = podfv. (5.16) 
We consider the following equation for Y 

— - kY" + K^Y = d^G in/x[0,T], (5.17) 
Po Po 

Y = on r X [0,r], (5.18) 
Y\t=o = pqUq in /. (5.19) 



where uq is given by 

5.4 Existence of a weak solution to the linear problem 
(I5.17p -( l5.19p by a Galerkin scheme. 

First we try to show that dfG is a function in £^(0, T; L^{I)). 



T fT r. 

6^112 _ / II q6/ ^ tPO \i\\\2 



/ P,^(-(^)')||i.(,) 

JO n n 



<C I Wpodfv" + dfv'Wl^d) + l.o.t 







(5.20) 



<C r \\pod!v\\HHi) + C f WdMrnd) + 1-o.t 
Jo Jo 

Now we begin our Galerkin scheme. Let {e„}„gN denote a Hilbert basis of 
Hi (/) . Such a choice of basis is indeed possible as we can take for instance the 
eigenfunctions of the Laplace operator on / with vanishing Dirichlet boundary 
conditions. We then define the Galerkin approximation at order n > 1 of (|5.19p 
as 1^ = Xir=o •^r(^)6ji ^iih A"(i) being the solution of the ODE system: 

Y o" 

{ — ,ek)mi) + {KY^,e'k)L^i) + (— efc)L2(j) = {dfG, ek)L^{i) 
^° ^° (5 21) 

Vfc e 0,...n 

Since each is in Hq{I), we have by high-order Hardy- type inequality p.ip 

that — G L^{I)- Therefore each integral in (15.2ip is well-defined. Further- 
Po 

more, as the {ci} are linearly independent, so are the { } and therefore the 

determinant of the matrix 



^/PQ y/PO J (ij)eN„ = {l,...,n} 
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is nonzero. This implies that our finite-dimensional Galerkin approximation 
(|5.2ip is a well-defined first-order differential system of order n+1, which there- 
fore has a solution on a time interval [0,T!„], where T„ may depend on the 
dimension n of the Galerkin approximation. 

Next we show that Tn > T, with T independent of n. 

Noticing that F„ is a linear combination of the ei{i G N„), we have that 

(^,r„)i2(,) -f + (^y„,r„)^.(,) = {dfG,Y„)L2(i) (5.22) 

Pq Po 
Hence, we have 

^ ^° (5.23) 



Y 

< PtG'||i2(/) -I- ||— =11^2(7) I |po||l°=(/) 



Using Poincare inequality ||Fn||^2(-/-| < II^t'IIl2(/) ^'^'^ Gronwell inequality, 
then we can find T > (independent of n) such that: 

1„ I ,,, , 1,2 



supC/^ + «/ ||r„||^^.(,) < / ||aJ>G||i2(,) + ||^||i2(,) (5.24) 
te[o,T] Ji Po Jo Jo VPo 

noticing (|5.20p and the fact v eCriM), 

sup C j^ + K f < A/-0 + CPdlTJlU^) (5.25) 

te[o,T] J I Po Jo 

where Ao is defined in (|5.8p . Thus, there exists a subsequence of (F„) which 
converges weakly to some Y G L^{0, T; Hq{I)), which satisfies 



sup C [ — + \\Y\\]j,,j)<No + CP{\\v\\x^) 

;e[o,Tl Ji Po Jo 



te[o,T] 

Now take the limit n ^ oo in (j5.2ip . we have 



(5.26) 



{-,ek)LHi) + n{Y', e't,)^!) + {^Y, ek)^!) = (5f G, efe)i2(,) (5.27) 
Po Po 

for every k. 

Hence, (15.19^ is satisfied in the sense of distributions, and that 

^* G L"(0,T;i/-i(/)) (5.28) 



Now we define 



Po 



ft 

Z= I Y{,T)dT + poU5, (5.29) 



W = [ Z{.,T)dT + poU4, 
Jo 



(5.30) 



and 

ft pti pt2 rts pt4 



x = f r r r r z{, T)dTdtidhdt2dti + v (5.31) 

Jo Jo Jo Jo Jo 

We then see that X e W^'^{[0, T];H^{I)) is a solution of ^J^-^J^, with 



dfx = Y. 



12 



5.5 Improving Space Regularity 

In order to prove that v G Ct{M) and then obtain a fixed point for the map 
Q :v 1-^ V, we need to estabhsh better space regularity for Z, and hence X and 

V. 

As Z defined in (I5.29p . then Z satisfies the following equation: 

— - kZ" + K^Z = d^G. (5.32) 
Po Po 

With high-order Hardy-type inequality, we have 
Zt Po 

Po Po (5.33) 

< \\y\\hhi) + IIpoIIl^||^||//m/) + WdtGWmi). 

So the regularity of Z = pod^v can be improved to L^{0, T; H'^{I)), and then 

V = — is well-defined and can be easily proved it is a solution to (I5.10p . 

Po 

Furthermore, as W defined in (j5.30p . we can see that W — podfv and W G 
L^{0,T; H^{I)). And we have a similarly estimate: 

Wt On 

>^\\W"\\mii) < W — Whhd + IRl^llffM/) + WdfGWmn) 

Po Po (5.34) 

< + ||p[,'||l~||W^||h^(/) + WdtGllnm). 

Hence that podfv € L^{0, T; H^il)) and dfv G L^{0, T; H^il)), and we have 

V e Xt- 

5.6 The existence of a fixed-point 

First it is clear that there is only one solution v £ L^{0, T; H^{I)) of (j5.10D with 
w(0) — uq, since if we denote by w another solution with the same regularity, 
then the difference 5v = v — oj satisfies Sv{-,0) — and poSvt — k[pqSv']' — 0, 
which implies 

a f PoSv^ + n f poSv'^ ^ (5.35) 



2dt Jj 



I 



which together with Sv(-,0) — implies Sv ^ 0. So the mapping v ^ v is well 
defined. 

Now we will prove v G Ct{M) when T is sufficiently small. 
First, we need to re-estimate L^(0, T; iJ^(/))-norni of Z = podfv. Like 
(j5.23p . we can easily have the following: 

and 

sup C f — + 2n r \\Y\\l,^j^< r \{d^G,Y)Lm)\ + \\^\\h(r) (5-37) 
te[o,T] J I Po Jo Jo \JPo 
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Since 



/ \{dtG,Y)L^i) 
Jo 




(5.38) 



<CTP{\\v\\%^) + CT sup 

tG[0,T] y/PO 

SO when T is sufficiently small, we can get 

sup C ( — + 2n I llFll^wn < AAo + CTP(||T7||2, ). (5.39) 
te[o,T] J/ Po Jo 

Considering (I5.32p . and using high-order Hardy- type inequality p.ip and 
the estimate (|5.39l) . we have 

c / II^"IIl2(j) 















Po 


Jo 


Y\\ 


HHI) + 






Jo 



(5.41) 



</ ll^lli..+ / (5.40) 
</ IIPolU~PII?fi(/)+A/'o + CTP(||T;||2.J 

<AAo + CrP(||F||^J. 
This implies 

\\Pod!v\\l^o,T;HHi)) < -^0 + CTPdl^JlP^J 

PMli=(o,T;ffM/)) ^^0 + CTP{\\v\\l^) 

The second inequality is following by using the high-order Hardy-type in- 
equality. The left part of Xt norm can be estimated in almost the same way. 
So finally we get 

M%^<Mo + CTP{\\v\\l^) (5.42) 

Take 

T < (5.43) 

- CP{M)' ^ ' 

we have \\v\\l;^ < 2Na- Let us fix M = 27Vo, then v G Ct{M). 

Now we have the mapping Q : v ~> v is actually from Ct{M) into itself 
for T — Tf^ satisfying (|5.43p . Then, we can get a sequence of functions S 
CriM), where — 9('y'^"^). It is obvious that i;*-"-* converges weakly in Xt- 

Furthermore, we have the following lemma which will show that pov^""^ converges 
strongly in P^(0, T; H^{I)) and hence w*^"' converges strongly in -L^(0, T; L^{I)), 
which will lead a fixed-point to system (|4.7p . 
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Lemma 5.1. For the sequence of functions u*^"^ we defined before, we 

have: 

(5.44) 

<CTF(||po(^^("+^)-^;("))|li.(o,T;/^M/))■ 
Proof. It is clear that — satisfies the equation: 

(„("+2)„„(n+l))|^^g=0. 

(5.45) 

Let U = po(w^"+^^ - similar as ([OS]) , we have: 

2 7/ Po 7/ Po 7/ 

< |([G(z;("+i))-G(«("))],{/)i.(,)| 

('n("+l)V2 _ ('n(")V2 

< C.ll /* po(«("+^^ - t'("^)'|li. + ^||po(t'("+') - t'("+^')'|li. 







<G^|| / po{v'^-+'^ ~v''-y)'\\l,+8C\\U\\l^. 





Then choose 5 small enough, and using Poincare inequality, Gronwell in- 
equality and high order Hardy type inequality, we finally have 

<CrP(||po(«("+^)-^^("^)|li.(o,T;/fi(/)))- 

□ 

Thereby, we prove the following Theorem: 

Theorem 5.2. // the initial data is smooth, then there exists a unique 
solution e Xt to the degenerate parabolic K-problem 7\ ) for sufficiently 
small T . 



6 Asymptotic estimates for v^^ independent of «: 

Our objective in this section is to show that the higher-order energy function E 
defined in (|2.25p satisfies the inequality 

sup E{t) < Mo + CTP{ sup E{t)) (6.1) 
te[o,T] te[o,T] 

where P denotes a polynomial function, and for T > taken sufhciently small, 
with Mq being a constant depending only on the initial data. The norms in E 
are for solutions to our degenerate parabolic K-problem (14. 7p . 
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According to (|5.2p . G Xt^ with the additional bound ||5^Uk||l2(q 7-^) < 00. 
As such, the energy function E is continuous with respect to t, and the inequahty 
(|6.ip would thus establish a time interval of existence and bound which are both 
independent of k. For the sake of notational convenience, we shall denote by 
V. We will generally follow the computation in [f?, sec. 6]. 



6.1 A /t-independent energy estimate on the S^^-problem 

Our starting point shall be the fifth time differentiated problem of (14. 5p for 
which we have, by naturally using dfv E L^{0,Ti^; H^{I)) (since v G Xt^) as a 
test function, the following identity: 



p,\ai-vf - I 3f iSiaf »' + « / pUsyf = 0. (6.2) 



^1 1-2 2^3 

Noticing the fact that d^v e L'^{0,T^; L'^{I)), which follows from (IPS)) 
and high-order Hardy-type inequality, (j6.2l) is well-defined. Upon integration in 
time, both the terms Ii and I3 provide sign-definite energy contributions, so 
we focus our attention on the nonlinear estimates required of the term X2 . 
We see that 



(6.3) 

Hence integrating (16. 2p from to t £ [OjT^], we find that 

~2 



Pod^v^iO) + / (dtvr^iO) - 3 / / {dfvrv'f^, (6.4) 
J/ Jo Ji V 

a=l ^ 

We next show that all of the error terms can be bounded by CiP(sup[Q E). 
First, it is clear that 

^dyfv'4<c fn^Wpodfv'wh 

Ji V Jo 

<C f \\v\\H2{\\po^tv\\],^+\\^tv\\l.) (6.5) 



< CtP{snpE). 

[0,t] 
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Then using integration-by-parts in time, we have that 
Jo Ji„^^ V 



J 



The term J can be written under the form of the sum of space-time integrals 
of the following types: 

Ji= j j Podtv'v'R{j/)podfv', 



J2 = J* l^pod!v'{v'fR{v')podtv', 
h = J^Pod!v'dtv'R{rj')podtv', 
Ji = J^pod?v'dtv'v'R{'n')podtv', 
J5 = j^Po{dyfR{v')podtv', 



Je = I* l^pod^v'{v'fR{v')podtv', 

J7= f I po{dtv'fR{v')podtv', 
Jo Ji 

Js^ f f po{dtvT{vTRW)Podtv'. 



IJ3I <C f \\pod!v'\\L.\\dtv'\\L.\\Rir,')\\L^\\podtv'\\L. 
Jo 



< CtP{supE), 
[o,t] 



and 



\J7\<C f \\dtv'\\U\R{v')\\L^\\podtv'\ 
Jo 



< CtP{supE), 

[0,t] 



(6.7) 



where R{ri') denotes a rational function of rj'. 
We first immediately see that 

<C [ \\v'\\Loo\\podtv'\\U\R{v')\\Loo < CtP{snpE). (6.8) 
Jo [0,t] 

Next, we have that 



(6.9) 



(6.10) 
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where we used Sobolev embedding inequalities in 1-D, || • 11^°= < QjH • and 

II • IIlp < Cpll • 11^1 , for all 1 < p < oo. 

J2, J4, J5, Jq and Jg can be estimated almost in the same way. 

1 * 
The term Jjb^d^—dfv'pldfv' 



can be estimated by 





E^{Mo + CtP{ sup E)) (6.11) 

t(E[0,T] 

in the similar way by using the fundamental theorem of calculus. 
Therefore, using Young's inequality, we have 



i / pod!v'{t)+ [ {dyf4(i)+K f I plid^Y < Mo+CtP{snpE), (6.12) 
^ Ji J I V Jo J I [o,t 



1 / pod^v^it) + I {podyfit) + ni I piidyf 



Ji 



and thus, using the fundamental theorem of calculus, 

/ pmv'f 

(6.13) 

< Mo + CtP (sup £;). 

[o,t] 

6.2 Elliptic and Hardy-type estimates for d^v{t) and v{t) 
Having obtained the energy estimate (|6.13p for the 9f -problem, we can begin 
our bootstrapping argument. We now consider the — i9?-problem of (|4.5p 



^-[d',^yjy,dy^ = -dy (e.u) 



which can be written as 

Po ?7 Po Po V Po r]'^ 

and finally be rewritten as the following identity: 



-2^[pldyni~^)~&pody^^. 

Po -n V 



(6.16) 



Here, ci and C2 are constants whose exact values are not important. 

Therefore, using Lemma p.2p and the fundamental theorem of calculus for 
the fourth term on the right-hand side of (j6.16p . we obtain that for any t S 
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[o,rj: 



2 

snp\\ — [pld^v']'\\L2 
[o.t] Po 

<sup||9MU.+sup||^[^°^]'|U. 
[o,t] [o,t] Po V 



+ sup||^[^]'IU.+sup||-[p2aM'IUHI3 f'^JlL^ 

\o.t] Po lo.t] Po Jo V 



(6.17) 



^2/3 

[o,t] " Po' V'^ [o,t] "Po Jo ^" 

77" 

We next estimate each term on the right hand side of (I6.17p . For the first term, 
we will use our estimate (j6.13D from which we infer for each t E [0, r„]: 

pllldfvf + \dtv'f]{t) < Mo + CtP{supE). (6.18) 

[o,t] 

Note that the first term of left-hand side of (j6.18l) comes from the first term of 
(j6.13p . together with the fact that dfv{t, x) — dfv{x, 0) + dfv{., x). Therefore, 
the Sobolev weighted embedding estimate p.2p provides us with the following 
estimate: 

/ \dtvf{t) < Mo + CtP{supE). (6.19) 

Jl [Q,t] 

The remaining terms will be estimated by simply using the definition of the 
energy function E. 

For the second term, we have that: 

II 1 rPodtVy_y.. 

< IKpoWIIlHI^IIl- + py[^]'IU^ 

<C||(po«0'llLHKlU- + lka^]llL= 

, II nPov".,, , .|| nPov'v".n 

+ H[^]\\L-+^\\Vt[^^]\\L^ 

<c\\{pou\y+ [\pov'uy\\L^\\v'\\hu+ f^wu ^^■'^^^ 

.In .In n 2 



+ C\\u', + / v[,\\l2\\v'\\IAK+ / v't 



+ C|K+ v[^\\L2\\pQv"\\''^^\\^UQ+ I ^Vt\\l2 
/*■ 1 

'II 2 



+ C||w |||fi||uo+ / ^^tll'ill / V ||i2||poUi+ / PoVttWm 
Jo Jo Jo 

< CsupEi{Mo + tP{supE)). 

[o,t] [o,t] 

Using the definition of E, then for any t € [0,Tk], we have 

supll — [^i^^]'||i2 <CsupE^Mo + tP{supE)). (6.21) 

[o.t] Po V [o,t] [o,t] 
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For the third term, we see that 

II r^o" i/M 
l| — [ — -J^l 
Po V 

^OII^'O^II , oil t^ PO^"'"' w , c;ll 'rW^^Viii 

+3||« +5||« [— 

<ciKiiiiii<+ /'^;iiiiii«o+ (6-22) 

Jo Jo " 

+ C'||pot'"||Ji||\/poWo + / VPow"|li2||uo + / 
Jo Jo 

+ C\\v'\\lA\u',+ f v'A\K\\u'o+ fvtWUW [' Pov"\\m, 
Jo Jo ^ Jo 

where we used the fact that || • ||^4 < Cp\\ ■ \\^^- Again, using the definition of 
E, the previous inequahty provides us for any t e [0,Tk] with 

supll — [^]'||l2 < Csup^i(Mo + iP(sup£;)). (6.23) 

[0,t] Po V [Q.t] [O.t] 

For the fourth term, we see that 



-[pldfv']'\\i^.\\3 [ < C[\\pody\\L. + \\dtv'\\i^.]tsup\\v\ 

Pa Jo V [o,t] 

< Ct P {sup E). 
[o,t] 



(6.24) 



Similarly, the fifth term yields the following estimate: 

IM,V^|U^(t)<C||po9*V|U^||r;"|U. 



<C\\po^fv'\\H^\\ I v"\\l- (6.25) 

< CtP{snpE). 

[o,t] 



Combining the estimates (I6.19p - (l6.25p . we obtain the inequality 

sup\\ — [pldtv']'\\L2 < Mo + CtP{snpE) + C sup E^Ma + tP{snpE)). (6.26) 

[Q.t] Po [0,t] [0,t] [0,t] 

We recall that the solution v to our parabolic K-problem is in Xt^ , so for any 
t e [0,T^],d'^v e H'^il). Notice that 

-ipld'^v']' = pod^ + 2p'odfv' (6.27) 
Po 

so (|6.26p is equivalent to 

sup\\pod^v" + 2p'od^v'\\L2 < CtP{supE) + CsnpEi{Mo + tP{snpE)) (6.28) 

[0,t] [O.t] [o.t] [0,f] 
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From this inequality, we would like to conclude that both ||9^w'||i2 and ||/Ooi9^^'"||l2 
are bounded by the right-hand side of (|6.28l) : the regularity provided by solu- 
tions of the K-problem allow us to arrive at this conclusion. 
By expanding the left-hand side of (j6.28p . we see that 



\\podtv" + 2p'^dtv'\\i. = Wtv"\\i.+A\\p'^dtv'\\i.+A J^podtv"p'odtv'. (6.29) 

We notice that the cross-term (|6.29p is an exact derivative with the regularity 
of dfv provide by our K-problem, 

Podfv"p'odfv' = 2l^pop'o§-\d?vT. (6.30) 

So that by integrating-by-parts, we find that 

podfv"p'odfv' = -2\\p'odfv'\\l, - J^podfv'p'^dfv', (6.31) 

and hence (j6.29p becomes 

\\pod?v" + 2p',d?v'\\l,^\\pod!v"\\h+2Ud?v'\\l,- I po5,Vpo9*V. (6.32) 



Since the energy function E contains \\,JpQd'^v{t)'\\[^2i^i'^, the fundamental the- 
orem of calculus show that 

Po9^Vo5^' < C\\^u'^ + / ^/p^d^^v'Wl. < Mo + CtP{snpE). (6.33) 

JQ [0,t] 

Combing this inequality with (j6.32p and (|6.26p . yields 



snv[\\pod^v"\\l^ + Udfv'\\h] 
[o,t] 

< Mo + CtP{su-pE) +Csui>E^Mo + tP{supE)). 

[0,t] [0,t] [0,t] 



3.34) 



and thus 



sup[||po9Mli. + \\p'od?v'\\l. + ||po9^'|li.] 

[o,t] 

< Mo + CtP{supE) + CsupEi{Mo + tP{supEj). 

[0,t] [0,t] [0,t] 

and hence with the physical vacuum conditions of po given by (|8.2p and 
we have that 

sup[||poa^"|li^ + l|5*V|li.] 

[0,t] 

< Mo + CtP{su-pE) + Csup£;3(A/o -I- tP{supE)). 

[0,t] [0,t] [0,t] 

which, together with (j6.19l) . provide us with the estimate 

sup[\\pody\\i. + \\dMl^] 

[0,t] 

< Mo + CtP{su-pE) + Csup£;3(A/o -I- tP{supE)). 

[0,t] [0,t] [0,t] 



3.35) 



3.36) 



3.37) 
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By studying the dxi — 9t)-problem of (14.51) in the same manner, we find that 
Po 

sup[||po^;"'|li. + lblll,.] 

[0,t] 

(6.38) 

< A/o + CtP{snY,E) + CsnpEi{Mo + tP{snpE)). 

[0,t] [O.t] [0,t] 

6.3 Elliptic and Hardy-type estimates for d^v{t) and dtv{t) 
We consider the — ^Qf-probleni of (14. 7p : 

^ [dt4]'-^ipldtv']'^-Vp-od!v, (6.39) 



Po ' ^ V'^' VPo 

By employing the fundamental theorem of calculus, it can be rewritten as 



Po VPo 

5 



w^+4=[^]'+^[^r (6.40) 



Po rj V 

for some constants ci and C2. 

For any t G [AjT^], Lemma p.2p provides the K-independent estimate 

2 

SUp||^[pg(9jS']'||L2 
[O.t] \JPO 



< sup 1 1 v^9M + sup 1 1 ^ [ ] ' I 

[OA [o,t] VPo 7?* 



+ sup||^[^2^]'|U.+sup||^[pga,VrilLH|3 /'^lU^ 
[O.t] y/Po [o,t] VPo Jo V 



(6.41) 



+ 6sup||p|a,V^IU.. 

We estimate each term on the right hand side of (|6.41l) . 

The first term on the right-hand side is bounded by Mq + CfP(supjo t] ^) 
due to (11331). 
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For the second term, we have that 



,4 J IlL^ 

Po ?7 



^ nVPodf^'ifov')' ^v'ipodfv'Y ^dfv'pQv'r] 

- " ^ " ^5 



<ciiv^4+/ ^d!v'\\L2\\{povy\\L 



+ C|k'||J,||u[,+ / v'AllAl / i;"|UHIPo'5tV||, 
■/o Jo 

<C\\^u'^ + f Va^9,V|U.||z;'|||,||<+ / u'"^ 



(6.42) 



ff7 



+ / ^d^v'\\L2\\pov"\\l,\\^u'a+ [ ^v['\\l, 

Jo Jo 

+ C\\vTm\Wo+ [' v',\\lA\y^u'^+ [\^d^v'\\i^2 



C\\v'\\hH+ v',\\1,\\poVp^u'^+ poVP^d!v"\\L2 



Jo 

+ C\\v'\\l,\\u', + f'v[\\l,\\ f v"\\L2\\podfv'\\H^. 
Jo Jo 



where we have again used fact that ||-||l°o<C||-|| a. Usmg the definition of 
E, it shows that for any t £ [0,T^], 



sup||^[^2^^]'|U. <Csup£;i(Afo + CiP(supi?)). (6.43) 
[o,t] V^'o ^ [OA [0,*] 



For the third term on the right-hand side of (I6.4ip . we have similarly that 

II 1 .Pldtv'^yu /,^ 



PO 11 

<2||(poWIU.|IVa^^IU-+5||Va^^IU2||po^"||lo 

rj'° ri'° 



<c\\{pou[y+ / ipodtvy\\L2\\y^jrodtvTi:2''\\iVFodtv'Y\\h-. 

Jo 

+ C\\d,v'\\U\ I'{pov")'\\l^ 
Jo 

<c\\{pou',)'+ f\podfvy\\L.\\u[+ [' d^v'wi-.^wiy^dtvy 



(6.44) 



■c\\d,v'r^.\\ I (poO'iu- 



3 — 3a 



where < a < i is given and < a < 1. 

The only term on the right-hand side of (|6.44p which is not directly contained 
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in the definition of E is \ \{^/po^tv'y\\22-a■ Then we notice tliat 

(6.45) 



1 



<\\^\\^2-.\\dtv'\\^.+\\^v','\ 



L2 



wfiere we have used the fact that 11 • ||lp < CII • 11 i , for all 1 < p < cx). So 

H 2 

(|6.44p and (|6.45l) provides us for any t e [0, r«] with 



sup||^[^2^]'||^, <C/supii;t(Afo + iP(supii;)). (6.46) 
[o,t] V/'o [o,t] [o,t] 



3 — 3a 

where < a = < 1. 

4 + 3a 

The fourth term on the right-hand side of (j6.4ip is easily treated as: 



Po - Jo V'* 

< C[\\pldfv"\\L2 + \\^d!v'\\L2]tsnp\\v\\H2 (6.47) 

[0,*] 

< CtP{supE). 

[o,t] 

Similarly, the fifth term is estimated as follows: 

\\pldfv'^\U.{t) < C||V^9,V|U.|I fpov"\\H^ 



V 

< CtP{supE). 
Combining the estimates (|6.42p - (|6.48p . we can show that 



(6.48) 



snp \\^[pld^v'y\\L2 < Mo + CtP{sup E) + C sup E'^ {Mo+tP{sup E)). (6.49) 

[0,t] y/PO [0,t] [0,t] [0,t] 

Now, since for any t g [0,Tk], solutions to our parabolic K-problem have the 
regularity d^v G H^{I), we integrate-by-parts: 

\\^[pld!v']'\\h = \\pld!v"\\l^+A\\y^op'od!v'\\h+2 f p'oPmvf]' 
y/Po Jl 



(6.50) 

Combined with (I6.49p . and the fact that podfv' = Jq Podfv' for the second 

term on the right-hand side of (|6.50p . we find that 

sup||p|9tV'||i2 <Mo + CtP{&npE) + CsnpE°'{Mo + tP{svipE)). (6.51) 

[0,t] [0,t] [0,t] [0,t] 



Now, since 



^ [pld^^v']' ^pld^y + 2^op'^dlv', (6.52) 
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the estimate (16.491) and (|6.5ip also imply that 
sup||yA)Po5tV||i2 <Mo + CtP{snpE) + CsupE°'{Mo+tPisupE)). (6.53) 

[0,t] [0,t] [0,t] [0,t] 



Therefore, 



[0,t] 

< Mq + CtP{supE)) + CsupE°'{Mo + tP{supE)) 

[0,t] [0,t] [0,t] 



(6.54) 



so that with (18.21) and 



snp[\\p-Sdfv"\\l^ + \\^d!v'\\U 
[o,t] 

< Ma + CtP{supE)) +CsupE"{Mo + tP{supE)) 

[0,t] [0,t] [0,t] 



3.55) 



Together with (j6.19p and the weighted embedding estimate p.2p . the above 
inequality shows that 



sup[\\pldfv"\\h + \\dMli] 



< Mo + CtP{supEj) + CsupE"{Mo + tP{supE)). 

[o.t] [o.t] [o,t] 



3.56) 



By studying the y/podx{ — 9?)-problem of (14. 5p in the same manner, we find 
Po 

that 



~ " ,."'l|2 I II a „,l|2 



sMMdtv'"\\i. + \\dtvl 

(6.57) 

< Mo + CtP{supE)) + CsupE"{Mo + tPisupE)). 

[O.t] [O.t] [0,t] 

7 Proof of Theorem 

7.1 Time of existence and bounds independent of k and 
existence of solutions to ( 12.221) 

Summing the inequality ((6l^ . (I07)) . ((OSl) . (|636l) . (|637)) . we find that 

snpE{t) < Mo + CtP{supE)) + CsupE°'{Mo+tP{supE)). (7.1) 

[0,t] [0,t] [0,t] [Q,t] 

As a < 1, by using Young's inequality and readjusting the constants, we obtain 

sup E{t) < Mo + CtP(sup E). (7.2) 

[o,t] [O.t] 

This provides us with a time of existence Ti independent of n and an estimate 
on (0,ri) independent of k of the type: 

sup E{t) < 2Mo. (7.3) 

[O.Ti] 

In particular, our sequence of solutions (vk.) satisiy the K-independent bound 
(|7.3|) on the ^-independent time interval (0,Ti). 
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7.2 The limit as k ^ 



By the k- independent estimate (|7.3I) . there exists a subsequence of which 
converges weakly to v in L2(0,r;7j2(/)). With r;(a;,t) = a; + J^v{x,s)ds, by 
standard compactness arguments, we see that a further subsequence of and 
77^ uniformly converges to v and 77', respectively, which shows w is the solution 
to t^rF3i - t^rM\i and v{x,0) = ua{x). 

8 The general case for 1 < 7 < 3 

If 7 7^ 2, we set loq — Pq^^, then physical vacuum condition shows that 

Wo > Cdist(x,a/), (8.1) 
when X € I near the vacuum boundary F, and 
du!Q 



> C when dix,dl) < a, (8.2) 
Wo > Ca > when d(x, £»/) > a. (8.3) 



Now we can use ujqv as intermediate variable and construct approximate 
solution to degenerate parabolic regularization just in a similar way to Section 
[SJ Noticing that the force term F ^ Jq Pody + C would not be smooth now, we 
need to require a certain high space regularity for it to keep the method described 
in Section [5] still work. Since we would require that X = luqv e H^{I), then 
from (jS.Sip and (15. ip . we will need the regularity that pody G H^{I). 

1 _A 1 

With (|8.ip and po = i^o , we will just require that ujq^^ G L'^{I)^ which 
means 1 < 7 < 3. 

So for 1 < 7 < 3, we can get the local wellposedness by doing the similar 
proof as 7 = 2 in Section [5] ~ [71 Details can be seen in [2] . 
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